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Abstract 



We report on a microscopic calculation of the current and magnetic field 
distribution in the mixed state of model s-wave and d-wave superconductors. 
For the d-wave case, we find that corrections to London theory are important 
at fields small compared to H c i- The field distribution is influenced by both 
superfluid-velocity dependence and non-locality in the current response func- 
tion of the mixed state. We compare our calculations with recent muon spin 
rotation measurements in high temperature superconductors. 
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The Meissner effect, in which weak external magnetic fields are exponentially screened 
from the bulk, is a basic property of superconducting metals. It is a consequence of the 
linear relationship between circulating currents (J) and the magnetic flux density (h) in a 
superconductor expressed by the London equation: 

h+^VxJ = 0. (1) 

In Eq.([l|) A, the penetration depth, determines the spatial extent of magnetic fields near 
the surface of a superconductor in the Meissner state. Measurements of this parameter in 
a particular material provide important information on the microscopic nature of its super- 
conducting state. For example, recent measurements of the temperature dependence of 
A, which microscopic theory relates to the quasiparticle energy spectrum 0, have provided 
important evidence for d-wave pairing 0] in high-temperature superconductors. 

Both nuclear magnetic resonance (NMR) and muon spin rotation (/xSR) measurements 
of A are based on the relationship between magnetic field inhomogeneity in the mixed state 
of a type-II superconductor and the penetration depth. In the mixed state, Eq. QT]) applies 
only outside the vortex cores and the London equation becomes, 

h + *E>L V xj = z$ ^ 5^ (r - R,), (2) 

c i 

where Rj is a vortex position and $o is a superconducting flux quantum. The right hand side 
of Eq.(||) accounts for the phase winding of the superconducting order parameter around 
each vortex. The magnetic field distribution inside the superconductor can be evaluated 
H from Eq.fl2|) and the Maxwell equation, V x h = (4n/c)J, once the vortex positions are 
specified. Eq.(^) should be valid provided that: i) the magnetic field is sufficiently weak 
compared to the upper critical field so that vortex cores occupy a small fraction of the 
volume; ii) the distance scales on which currents and fields vary inside the superconductor 
are larger than the non-locality length of the current response kernel; iii) and the fields 
and currents are small enough to be in the linear response regime. This last condition is 
especially pertinent in the case of d-wave superconductors since nodes in the energy gap 
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lead to large non-linearities ||[7| in the current response, even at low temperatures. At 
present there is no detailed theory which accounts for non-linear supercurrent response 
in the inhomgeneous mixed state. In this article we report on exact microscopic mean- 
field-theory calculations of the current and magnetic field distributions in the mixed state 
of model s-wave and d-wave superconductors. We summarize our results in terms of an 
effective penetration depth which would be inferred from the calculated field distribution 
under the assumption that Eq.(@) applies locally. We find that this effective penetration 
depth increases with increasing field, and hence increasing supercurrent densities, much 
more strongly for d-wave superconductors than for s-wave superconductors and that the 
linear temperature dependence of A which signals <i-wave superconductivity in the absence 
of a magnetic field saturates at a temperature k B T* ~ (H/ ^ c2 ) 1//2 A . 

Our numerical calculations start from self-consistent solutions H-jlTf of the mean-field 



Bogoliubov-de Gennes (BdG) equations flT2| , |T3l1 for generalized Hubbard models on stacked 



magnetically coupled square lattices. We regard the Hubbard models as plausible low-energy 
effective models for high-temperature superconductors. The class of Hubbard models we 
consider is paramaterized by the nearest neighbor hopping energy (t), the on site interac- 
tion energy (U) and the nearest neighbor interaction energy (V) and the parameters can be 
chosen so as to obtain s-wave or d-wave superconductors at zero magnetic field. The model 
parameters for the calculations discussed below were chosen to give low temperature coher- 
ence lengths comparable to estimates made for typical high temperature superconductors. 



We have solved the BdG equations [14] self-consistently for s-wave and d-wave mixed states, 



assumming [15] a constant magnetic flux density Hz, as detailed in an earlier publication. 



16| The current flowing along the bond of the lattice linking site i and site j is related to 



the self-consistent BdG quasiparticle amplitudes by 

Jij = ^ im {^ n«r«i/« + v >?{i - /„)]}, (3) 

where a is in-plane lattice constant, d is the separation between square lattice layers, and 
f a = [exp(£' Q //csT) + 1] _1 is the Fermi function for quasiparticle excitations of energy E a . 



The spatial variation in the flux density h in each plaquette of the square lattice is then 
evaluated from a discrete version of the magnetostatic Maxwell equation: 

, . , 47rz „ G x J(G) , , AS 

h(r) = rT^ ex P( zG ■ r ) + H - ( 4 ) 

C G^O U 

Typical results are shown in Fig. [I] for both s-wave (left column) and d-wave (right 
column) superconductors. In these figures we show the spatial dependence of the current 
density, the local magnetic flux density, and a local penetration depth which we define by 

2 _ 4tt Vx J 

These results were obtained for T = at fields with approximately 900 plaquettes per 
superconducting flux quantum which corresponds, for the model studied, to H ~ 0.04H c2 
in both s-wave and <i-wave cases. 

The screening current circulating around the vortex core is largest in magnitude at the 
core edge. Qualitatively, this property can be easily understood in a Ginzburg-Landau type 
of analysis, where the current density is related to the order parameter ip — \tp \ exp(i<f)) and 
the superfluid velocity v s by 

J = 2e|Vfv s = ^H 2 (W0--A). (6) 

m* v c ' 

Since (f) winds by 2ir around a vortex and \if)\ vanishes linearly at vortex center, J ~ r near 
the core center and J ~ r _1 (for r < A) outside the core. Our results show a distortion in 



the circular flow pattern generally expected around the vortex core which we attribute [17 



to a tendency for current to flow along the principle symmetry axes of the square lattice. 



We find that the local penetration depth is tolerably constant outside the vortex core [JTS 
as assummed in the simple London theory, except where the current flow bends relative 
abruptly. Similar band structure effect has been seen in scanning-tunneling-microscopy 
experiments for conventional superconductors by Hess et al. |I9" |. However, near and inside 
the vortex core Xj^, changes rapidly, presumably principally because the magnitude of the 
order parameter is changing rapidly. There is no qualitatively difference between s-wave and 



ci-wave cases apparent in the results at a particular temperature and field shown in Fig. [l|. 
The deep minimum in the local penetration depth in the vortex core is effectively replaced 
by a 5-function in the London approximation. 

In the Meissner state of a rf-wave superconductor, the dependence of the penetration 
depth on superffuid velocity leads to important non-linear terms in the current response and 
to a field dependent effective penetration depth. In the mixed state, the dependence of the 
magnetic field on two spatial coordinates and the presence of vortex cores where the order 
parameter magnitude is not constant complicate the analysis of non- linear effects. In fact, 
the current density in the mixed state is largest near the vortex cores and this complicated 
region might therefore be expected to contribute importantly to the field dependence of any 
effective penetration depth, making the development of an analytic theory truly difficult. 
The magnetic field at any point in the vortex lattice depends on the current distribution 
everywhere and hence its distribution function depends on some complicated average of the 
non-local, non-linear, and inhomogenity effects which appear in the local penetration depth. 
It is useful to define a magnetic-field-dependent effective penetration depth A in the mixed 
state by mimicing the analyses of /iSR or NMR experiments which are based on the London 
equation. This A may be thought of as an appropriate average of the local penetration depth 
discussed above. 

In a /iSR experiment one measures the time evolution of the muon polarization vector 
precessing about the internal magnetic field. The frequency spectrum is proportional to the 
magnetic field probability distribution in the sample. When the London equation is valid 
nearly everywhere it can be shown [2(| that A is related to the width of the magnetic field 
probability distribution 

P(h) = jjdr6(h-h(r)) (7) 

by the following expression, 

A" 4 = C (Ah 2 ) = C Jdh P(h) (h - H)\ (8) 
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where the constant C is dependent on the structure of the lattice. We define the effective 
penetration depth of a vortex lattice by this equation using the exact field distribution. This 
definition permits a direct comparison between theory and experiment. 

In Fig. |2| we show the effective penetration depth A e //, as a function of temperature. We 
find that A e // increases with increasing magnetic field for both d-wave and s-wave supercon- 
ductors, but much strongly in the <i-wave case. Even in a relatively weak field, H ~ 0.04if C 2, 
the low-temperature behavior in two cases is very different. In the s-wave case, A e // essen- 
tially preserves an activated T-dependence. On the other hand, the linear temperature 
behavior for the (i-wave A in the Meissner state saturates at T*(H), and crosses over to a T 2 
dependence. The crossover temperature, T*, can be crudely estimated by noting that the 
typical local pairing momentum in the vortex lattice state is ~ h/£ [plj] , where 27r£ 2 H = <3>o- 



This local superfluid velocity changes the quasiparticle density of states ||[7,21,22] and hence 



A for energies below kBT*(H) ~ vpfoji ~ A ^ H/H c2} where A is the maximaum quasi- 
partice energy gap in the Meissner state. This estimate of the crossover temperature is 
consistent with the numerical results shown in Fig. ^| and similar results obtained at still 
stronger magnetic fields. 

Our calculation should be compared with recent /zSR experimental data on YBa2Cus06.95 
reported by Sonier et al. 0. These authors find that the penetration depth increases with 
field. The increase of approximately 15% they find at the strongest field studied, H = 
2 Tesla, should be compared with the approximately 30% increase we find for H/H c2 ~ 
0.04. (H c2 is belived to be ~ 200 Tesla for YB^CuaOg.gs). This is consistent with the 
expected approximate relationship between the zero-temperature penetration depth and 



the quasiparticle density of states at the Fermi energy which is oc yH ||lq , |21] , |23l| . The 
observed change in X e ff at this field is much larger than would be expected for an s-wave 
superconductor, adding to the evidence for (i-wave pairing in cuprate superconductors. This 
experimental study has insufficient data at sufficient low temperature and sufficient strong 
fields to verify the crossover T 2 behavior predicted here. 

This work was supported by the National Science Foundation under grant DMR-9416906. 
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FIGURES 

FIG. 1. Top row: Current circulating around a vortex. Arrows represent the directions of flow. 
The length of an arrow represents the magnitude of the current at the site. Middle row: Local 
penetration depth as defined in the text. The minimum in this quantity occurs at the center of 
the vortex core. Bottom row: Magnetic field distribution. The peak occurs at the center of the 
vortex core. The plotted quantity is the difference between total field and the spatial average field 
in arbitrary units. The applied field is H = 0.0AH C 2. The figures in left and right columns are for 
s-wave case and d-wave case respectively. 

FIG. 2. Temperature dependence of inverse square of the effective penetration depth X e ff(T) 
in the mixed state (H ~ O.OAH^) for ti-wave and s-wave superconductors, normalized by the 
zero-field, zero-temperature penetration depth Ao(0), is plotted against reduced temperature T/T c . 
H n represents the field corresponding to the inter-vortex spacing of na. For typical Cu-Cu distance 
in high T c materials, H 30 and H 2 g, are 8.0 and 9.2 Tesla, respectively. The zero-field penetration 
depth calculated with the same models are shown for comparison. 
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FIG.l 
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